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An accurate knowledge of the coalescing binary gravitational waveform is crucial for experimental 
searches as the ones performed by the LIGO- Virgo collaboration. Following an earlier paper by the 
same authors we refine the construction of analytical phenomenological waveforms describing the 
signal sourced by generically spinning binary systems. The gap between the initial inspiral part 
of the waveform, described by spin- Taylor approximants, and its final ring-down part, described 
by damped exponentials, is bridged by a phenomenological phase calibrated by comparison with 
the dominant spherical harmonic mode of a set of waveforms including both numerical and phe- 
nomenological waveforms of different type. All waveforms considered describe equal mass systems. 
The Advanced LIGO noise-weighted overlap integral between the numerical and phenomenological 
waveforms presented here ranges between 0.95 and 0.99 for a wide span of mass values. 
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I. INTRODUCTION 

The experimental program for gravitational wave detection is on the way, as the network of kilometer-scale interfer- 
ometers formed by the Laser Interferometer GravitationaFwave Observatory (LIGO), Virgo and GEO are presently 
in science run or undergoing substantial upgrades taking them to their advanced version [IHS] . 

Among the possible sources one of the most promising is represented by the coalescence of compact binary systems 
of neutron stars or black holes. The coalescence of binary systems is usually described in terms of three distinct phases: 
the inspiral, the merger and the ring-down. The inspiral phase allows for an accurate analytical description via the 
so-called Post-Newtonian (PN) expansion, see for instance [3] for a review. The ring-down also admits a perturbative 
analytical model, as it describes the damped oscillations of the single object resulting from the binary coalescence, 
as a superposition of black- hole quasi- normal modes [5j . The merger phase is however fully non-perturbative and for 
generic systems it has not been described analytically but rather by numerical simulations. 

During the last six years numerical relativity has made tremendous progress in describing the full coalescence of 
a binary system beginning with [BHH] and more recently [51115) . see ,16nr8] for reviews, and it can now produce 
waveforms for generic spin orientations, with moderate spin magnitude (< 0.9) and mass ratios (< 10 : 1). 

Match-filtering techniques represent a powerful tool for seeking signals, but they need a detailed knowledge of the 
waveform in order to have high efficiency and to be useful for parameter estimation. They are extensively used in 
LIGO- Virgo experimental searches in order to uncover weak signals buried into noise, see e.g. [T5] for the results of a 
recent match-filtered search. In such searches real data are compared with banks of template waveforms, made of a 
large number (tens of thousands) of templates: due to the computational cost of numerical simulations, it would be 
impractical to numerically generate the waveforms necessary to populate template banks. 

There exist nowadays analytically constructed waveforms describing the entire coalescence of a binary system. 
They have been achieved in the Effective One Body (EOB) construction [201123] for non-spinning systems, and in the 
EOB-spin waveforms [24] . for binaries with spins aligned with the orbital angular momentum (thus non-precessing 
spins). In both cases of spinning and non spinning waveforms, a comparison with numerically generated waveforms 
is needed by the EOB construction in order to calibrate some free parameters of the model. Another method 
for generating analytical waveforms for spinning non-precessing binaries is by joining PN-generated inspiral with 
a fit of numerical waveforms to construct ■phenomenological waveforms, as done in |13j . In this paper, we further 
investigate on the family of phenomenological waveforms introduced in [25!|, dubbed PhenSpin, which are analytical 
waveforms describing the entire coalescence of generically spinning compact binaries. In particular a restricted set of 
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phenomenological waveforms [13] describing non-precessing systems are used here together with the same set of fully 
precessing numerical waveforms used in |25j to tune this new version of PhenSpin waveforms. 

The PhenSpin waveforms have been constructed by joining the perturbative PN description of the inspiral to the 
ring-down phase by a phenomenological phase which plausibly describes the evolution of the waveform in between. 
With respect to the previous work introducing these waveforms, here we give a slightly modified (improved) version, 
identical in spirit but better tuned in some technical aspects, and produce new results assessing their faithfulness 
to numerical simulation in presence of detector noise. The improved details allow to obtain a slightly better match 
between the PhenSpin waveforms and the set of test waveforms they have been compared with. 

Waveforms describing generically spinning coalescing binaries are not suitable for searches employing match-filtering, 
as the size of a template bank increases exponentially with the number of template parameters: since spinning 
waveforms depend on several parameters (masses, spin components of the binary constituents, angles defining the 
orientation of the source with respect to the observer) it is not practical to construct one template bank to cover the 
entire spin parameter space. Non spinning, or at least non-precessing waveforms, are usually preferred for template 
bank construction, with the exception of the so-called Physical Template Family, representing a single spin family 
waveform V2& , which however can effectively describe also doubly spinning physical systems [271 . 

The availability of generically spinning waveforms is however badly needed to assess the efficiency of experimental 
searches based on banks of non-spinning templates, like |19j . Moreover fully spinning waveforms can be used as 
templates in connection with parameter estimation via Bayesian inference methods, which can be used as follow-up 
analysis to perform searches in the parameter space with full dimensionality, but restricted to a small subset of the 
entire space, as determined by lower level triggers. 

The paper is organized similarly to [25] and the exposition has been kept here as self-contained as possible. In sec.[Tl| 
the analytical waveform construction and the waveforms used for calibration are revisited. Differences with respect 
to the old version of PhenSpin are described. In sec. |III| the results are presented, in the form of comparison between 
analytically and numerically generated waveforms, which reproduce the dominant quadrupolar mode I = m = 2. In 
sec. 



IV the conclusions that can be drawn from the present work are reported. 



II. THE METHOD 



Following the original introduction of PhenSpin waveforms given in |25j , the present work revisits the construction 
of analytical gravitational waveforms generated by the coalescence of spinning binary systems. The waveforms used to 
construct and calibrate our analytical model include the numerical waveforms used in the previous paper, describing 
equal mass binary systems (m = mi = 7712), with spin magnitudes |Si| = IS2I = O.dm? and starting with S2 
orthogonal to the initial orbital angular momentum (where Si_2 denote the binary constituent spin vectors and we 
posit Gn = c — 1). Here in addition, to cover a portion of the parameter space not addressed by the numerical 
simulations, we use also four phenomenological waveforms of the type described in |13| . generated via the LAL 
libraries |29j to calibrate the PhenSpin in the case of aligned spins, which together with the numerical relativity 
waveforms form our set of test waveforms. 

The description of the dynamics adopted here models the inspiral phase via the standard TaylorT4 PN formulae, 
see [30' for definition and comparison of different PN approximants in the spin-less case. Un-like the non-precessing 
case, the knowledge of the time- varying amplitude and phase is necessary but not sufficient to determine the waveform 
from spinning precessing binaries, as it must be complemented by the spin and angular momentum evolution, see e.g. 

£ = -^(s, + S2), 

where 

^^1,2 ^ (4 + i T ^<5) Ln (2) 

and v = mim2/M^ is the symmetric mass ratio and S = (mi — 1712) /M, being M = mi -f m2 the total mass of the 
binary system and mathbfL the orbital angular momentum unit vector. It is convenient to define an orbital phase 
4> = J cJorb dt whose evolution is given by 



(3) 
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where F{v) and E{v) are respectively the flux emitted and the energy of a circular orbit with angular frequency tOorb, 
related to the main gravitational wave frequency few via few = ^orb/T^- 
By parametrizing the orbital angular momentum unit vector L as 

L — (sin L cos a, sin t sin a, cos t) (4) 

it is convenient to introduce the carrier phase ^' given by 

da 

— =a.„.,-cos.-. (5) 

Numerically generated waveforms are usually decomposed in spherical harmonics, in particular the five quadrupolar 
modes {I = 2) are the only non-vanishing at the lowest order in w, and the I — 2,m — ±2 mode are the dominant 
ones. As determined by the PN analysis, the I = 2,m = 2 mode in the inspiral phase i^^^ '^^^Y one which 

will be used here for comparison with test waveforms, which can be expressed in terms of the usual plus and cross 
polarizations as ft.2*2*^^ — ^ ihx) is given by formulae which can be found e.g. in [31j . which here we re- write in 
the following form 



f^r'it) = -2/^^-' [(1 - civf^m - 55.)) (cos^ (./2) + sm^(./2)e^^(*-")) + 



S 

V- sm 1 1 sm"" -e 
3 



fsin2-e'(*-2")+cos2-e'(*+2")') 
V 2 2 J 



(6) 



+ 0(«5), 



where spin-dependent terms in the amplitude have been neglected as well as terms of order higher than w^, t- 
dependence is understood in 5* and a and the z— axis used for defining Z,m modes is parallel to the initial total 
angular momentum. 

Note that differently from |31) and the previous version of the PhenSpin waveforms, the v'^ terms have been added 
weighting them by a phenomenological function c{v) defined as 



c{v) 



exp[-(l-0.05/a;orf,)V2] Wo^b < 0.05 

1 UJorb > 0.05 



which has the role of turning corrections on at values of the orbital frequency AIu]orb ?J 0.05, as otherwise poor 
matching with numerical simulations would be obtained. 

The m = —2 mode can be obtained via ft,2,-2(^) = ^2 + ^^id in the equal mass case /i2,-2 = ^22 holds: we 
thus focus the calibration of our phenomenological model on the /i2,2 mode. 

The functions F{v) and E{v) are necessary to determine the orbital phase and they are known up to 3.5PN order 
as far as orbital effects are concerned, and up to 3PN and 2PN level for respectively S1.2L and S1S2, SiSi, S2S2 
interactions, see [32-34J for recent derivations of spin-orbit and spin-spin interaction effects. 

According to studies in the non-spinning case [21] [35l [36] , the Taylor T4 appears to be a very good approximant up 
to a frequency u) — it few — 0.1/M for the equal mass case, even though its faithfulness seems to worsen for different 
mass-ratios (which however are not considered in this work). |45| 

The PN evolution ([6]) is halted a.t t = tm, when uiorb reaches the value Wm that is determined by comparison with 
the test waveforms. For Worh > '^m i'-^ort is monotonically increasing) the angular frequency is evolved according to 

^orb{t) — Z h LOq, UJm < UJorb < ZUJrd and t,n < t < t^ac , (7) 

l-t/lA 

where the three unknown parameters uq^i and Ta are determined by requiring continuity of coorb and its first and 
second derivatives at the matching point defined by cjrn- 

The damped exponentials describing the ring-down phase are attached at the instant of time trac when ujorb reaches a 
fraction z of the ujrd value, the specific value of z has been determined like by comparison with the numerical rela- 
tivity waveforms, as described below, and differently with respect to the previous version of the PhenSpin waveforms, 
where z has not been fit to numerical simulations but kept constant to a convenient value. 

Differently from the first PhenSpin version [23], where it has been kept constant at their value at t — tm, the 
angular variable a is evolved with a similar phenomenological formula 

da ai 

■ "0, (8) 



dt 1 - t/Z 



A 
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where the parameters di.o are determined by requiring the continuity of a up to its second derivative, and Ta is the 
same as determined in eq.Q. 

Finally the usual ring-down description of the waveform is used 

^2'2^(i) - E^"'^'"^"^'"'"'"* *>W, (9) 

n 

where it is used that the ring-down phase is described by adding damped exponentials of increasing inverse damping 
time, the overtones, with complex constants ^„'s. Here we assume that given the moderate spin values we are 
considering, the direction of the final spin of the black hole is parallel to the initial total angular momentum. We 
have checked that during the PN-inspiral phase (t < t^) this is indeed the case to very accurate precision (better 
than 10~*) for all spin configurations considered here. In the same numerical simulations used here are analyzed 
and a maximum misalignment angle 9 between the final spin and the initial total angular momentum is found to be 
around 9 ~ 0.24rad ~ 13° (see fig. 3 of |42]). Further investigations are necessary to assess the importance of this 
effect on the actual waveform shape. 

Allowing more overtones requires to fix more coefficients, which can be done by admitting continuity of the waveform 
to the appropriate level: using n overtones requires matching the waveform up to its 2 x (n— l)-th derivative, as each 
overtone involve the determination of a complex (or two real) constant parameter(s). The construction that inspired 
our work is the FOB matched to numerical relativity waveforms (usually referred to as EOBNR), introduced in [5T] 
where the waveform is assumed circularly polarized (i.e. /ix(5' + 7''/4) = ±/i-|_(^)), so that the real and imaginary 
part of An for the ^ = 2, m = 2 are not independent parameters. Here however we do not assume circular polarization, 
as in general terms of order in eq.([6| (for unequal mass systems) will spoil this property and the stitching of the 
ring-down modes is performed independently on the real and imaginary part of each multipolar mode. For any such 
mode defined by a (l,m) pair there is an infinity of overtones with increasing damping factors, but for our practical 
purposes retaining only two overtones is enough. 

As described in |43j each overtone with given /,to will be in general a superposition of the two modes which are 
usually designated by I, m and I, —m. Here we stick to the prescription adopted in [3T] where only the m > mode is 

stitched to the inspiral waveform. We have verified that adding the I — 2,m ~ ~2 mode to h'^!^' would not improve 
the fit to the numerical waveforms, at the expense of introducing additional constant parameters which have to be 
fixed by imposing further continuity requirements. 

The values of the ring-down frequencies and damping factors of the three lowest overtones of the I < 4 modes can 
be read from [37] as a function of the mass and spin of the final object created by the merger of the binary system. 
We estimate the final mass by taking the algebraic sum of the constituents' masses and the negative binding energy 
once ujrn is reached, and the final spin according to the phenomenological formula given in eq.(5) of [38J. 

The analytical waveforms just described have been quantitatively confronted with the set of test ones by computing 
the noise-less overlap integral 

h,M^'^ {hi{f)m) + K{f)h^{f)) df (10) 
maximized over initial phase and time of arrival, where normalized waveforms have been considered 

The angular frequency a;^ and the parameter z have been determined by comparison with a first set of short test 
waveforms (4-6 cycles long) by picking the values maximizing the overlap integral ( |10[ ) (with a precision respectively of 
±5 • 10~^/M and 0.01). The set of short waveforms have initial orbital frequency ujorb-ini ~ 0.05/Af (for comparison, 
<^rd ^ 0.3/M). Note that despite ujorb;ini being not too far from the values of oJm obtained through the fit, see tab.jlj 
it still allows at least one oscillation cycle before the onset of the phenomenological phase for all test waveforms. 

The numerical waveforms in the set of test waveforms have been generated with MayaKranc. The grid structure 
for each run consisted of 10 levels of refinement provided by CARPET [391, a mesh refinement package for CACTUS 
[40] . Sixth-order spatial finite differencing was used with the BSSN equations implemented with Kranc |41J. The 
outer boundaries are located at 317Af and the finest resolution is M/77. Waveforms were extracted at 75M. A few 
waveforms were generated at resolutions of {M/64, M/77, M/90}, and convergence consistent with our fourth order 
code is found. The short (long) runs showed a phase error on the order of 5 • 10~^ (5 • 10~^) radians and an amplitude 
error of « 2% (« 5%). 

The numerical waveforms consist of two sets: the first set consisted in 24 few-cycle-long waveforms, representing 
mostly the merger and ring-down phases of a coalescence which, together with four phenomenological ones with 
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aligned spins, has been used as described above to fix the values of and z for the corresponding values of initial 
spins. All of the numerical waveforms have initial spin ^^I'tt? = (—0.6,0,0) in the reference frame in which the 
initial L ~ (0, 0, 1). The different values of the first dimension-less spin have been obtained by rotating the (0, 0, 0.6) 
vector by 15 degrees in the x-z plane. This set of numerical waveforms has been completed by the addition of four 
phenomenological waveforms: one spin-less and three with spins aligned with the orbital angular momentum L and 
the same magnitude as above (one with both spins aligned with L, one with spins pointing in opposite directions and 
a third with both spins anti-aligned with L). Once determined the values of ujm and z for each waveform, their values 
for generic spins have been determined by assuming an analytical dependence on the dimension-less spin xi.2 defined 
as xi,2 = Si^2/'7ii 2i according to 

MuJrn = ao + ai(xi^ + X2z) + a25{xiz - X2z) + a3(xiX2) + 

a-iixl + xl) + a^{xizX2z) + a^iXiz + xL) + {xiz + X2^)x (11) 
[a7(xiX2) + asixl + xl) + a'i{xizX2z) + a-wixlz + xD] ■ • ■ , 

where the spin components are understood in a frame where the orbital angular momentum is along the z axis and 



higher powers of the spin components have been neglected. From eq. (11) one can note that since spins are evolving 
in time, LOm is also slightly changing with time, the explicit values taken by the coefficients are reported in tab. |ll] 
The Xi values used to fit are taken at t = t^. 

The dependence of the Worb evolution equation on L and Si ^2 implies that the spin components parallel to the 
orbital angular momentum enter already at linear level, whereas the dependence on the spin components in the plane 
of the orbit starts from the quadratic level. The coefficients may depend on the symmetric mass ratio but it is 
assumed here that they can be analytically expanded around their value at (5 = 0, according to 

a,{5) ^ a, + S 4^^^ + 6\l^'> + . . . . (12) 

Note that anti-symmetric combinations of spin components do not appear for (5 = 0. Given the specifics of the test 
waveform set we used (all having S ^ 0) we could not calculate the coefficients a^'^ nor 02. 
An analog formula has been assumed for z 

bn + bi{xiz+X2z)+b25{xiz-X2z) + b3{XiX2)+ ^^^^ 
hixl + xi) + b5ixizX2z) + beixlz + xlJ , 

where terms cubic in the spins have not been necessary here. Results from the fit of z are reported in tab. 

These values have then been tested by computing the faithfulness of the now fully calibrated PhenSpin wave- 
forms with respect to a second set of long waveforms, consisting of 8 long numerical waveforms (12-15 cy- 
cles long) plus the 10 phenomenological ones: 4 with the same parameters as above (initial oJorb — 0.03/M) 
and six additional ones characterized by spins aligned with the angular momentum and chi-pair values give by 
(xi^,X2^) = {(0.3, 0.3), (0.3, -0.3), (-0.3, -0.3), (0.8, 0.8), (0.8, -0.8), (-0.8, -0.8)}. The faithfulness of a pair of 



waveforms (/ii, /12) is quantified by the noise- weighted version of the overlap integral 10 which we rewrite in terms of 
e, the mismatch parameter (see e.g. [44] ) 



hiifmf) + hi(f)h2{f) 

Snif) 



Snif) J \J SM) 



df] / ^^^df 



(14) 



where maximization over initial phase and arrival time is understood and Snif) is the single-sided power spectral 
density of Advanced LIGO strain noise and waveforms are normalized with respect to noise weighted integrals. 
The results of the above integrals comparing PhenSpin waveforms and long test waveforms with identical physical 
parameters are described in the next section. 



III. RESULTS 



The analytical waveforms have been calibrated by comparison with 28 short test waveforms described in the 
previous section and the results obtained by maximizing the overlap given by eq. ( 10 1 are reported in tab. ffl The 



6 



determination of the LOm's giving the best overlap for different spin values allowed to evaluate some of the coefficients 
in the phenomenological formulae (11) and (13), as given in tab. In) 

Once fixed the value of and z for generic spin values, it is possible to generate analytical waveforms with any 
specific initial condition without any further tuning: the value ofuim,z will be determined analytically via eqs. (11) 
and ( 13 1 with the unknown coefficients arbitrarily set to zero. It is then possible to generate waveforms with no 
tunable parameters for comparison with the second set of long waveforms. The results of the faithfulness integrals 
described in eq. ( [l4| are reported in fig. [ijfor the range of masses 60-100 Mq. The lower range corresponds to the 
minimal mass value enabling the long waveforms to start at a physical frequency which is smaller than the lower edge 
of the sensitive band (which we assume to be around 20Hz) whereas beyond the upper range value only the ring-down 
phase is in-band (for reasonable distances of the sources). 



IV. CONCLUSIONS 



We presented an analytical method to produce complete gravitational waveforms from spinning coalescing binaries. 
The free parameters of the model are the values of the orbital frequency at the transition from the inspiral to the 
phenomenological phase and at the transition from the phenomenological to the ring-down phase. After a calibra- 
tion process involving the dominant multipolar mode obtained by numerical relativity and other phenomenological 
construction from a different family, all the parameters have been fixed and the PhenSpin are ready to use, once fed 
with the physical parameters of the coalescing binaries (mass, spins, inclination angles, initial phase). We computed 
noise-weighted overlap integral obtaining values between 0.95 and 0.98 for a wide range of masses. Further investiga- 
tions are necessary in order to assess how such mismatch can affect the precision of parameter estimation performed 
via these family of waveforms. 



Acknowledgments 

It is a pleasure to thank the organizers of the NRDA conference in Waterloo for the stimulating scientific environment 
created by the meeting. This work was supported by the FAI grant from INFN and by the NSF grant PHY-0653550, 
PHY0955773, PHY-0925345, PHY-0941417 and PHY-0903973, PHY0955825, TG-PHY060013N. 



References 



B. Abbott et al. [ LIGO Scientific Collaboration ], Kept. Prog. Phys. 72 (2009) 076901. |arXiv:0711.304l] [gr-qc]]. 
F. Acernese, M. Alshourbagy, P. Amico, F. Antonucci, S. Aoudia, K. G. Arun, P. Astone, S. Avino et al., Class. Quant. 
Grav. 25 (2008) 184001. 

B. Willke [ LIGO Scientific Collaboration ], Glass. Quant. Grav. 24 (2007) S389-S397. 
L. Blanchet, Living Rev. Rel. 9 (2006) 4. 
S. A. Teukolsky, Astrophys. J. 185 (1973) 635. 

F. Pretorius, Phys. Rev. Lett. 95 (2005) 121101 [arXiv: "gr-qc/0507014". 

M. Campanelli, CO. Lousto, P. Marronetti, and Y. Zlochower, Phys. Rev. Lett. 96 (2006) 111101 [arXiv:gr-qc/0 511048 |. 
J.G. Baker, J. Centrella, D. Choi, M. Koppitz and J. van Meter, Phys. Rev. Lett. 96 (2006) 111102 arXiv:gr-qc/0511103'. 
F. Herrmann, I. Hinder, D.M. Shoemaker, P. Laguna and R.A. Matzner, Phys. Rev. D 76 (2007) 084032 
[ arXiv:0706.2541 ] . 

W. Tichy and P. Marronetti, Phys. Rev. D 78 (2008) 081501 [arXiv:0807.2985 [gr-qc]]. 

M.A. Scheel, M. Boyle, T. Chu, L.E. Kidder, K.D. Matthews, H.P Pfeiffer Phys. Rev. D 79 (2009) 024003 [arXiv:0810.1767| 

[gr-qc]]. 

[12] J. G. Bake r, W. D . Boggs, J. Centrella, B. J. Kelly, S.T. McWiUiams and J.R. van Meter, Phys. Rev. D 78 (2008) 044046 
rarXiv:080 5TT428] [gr-qc]]. 



7 



# 


Overlap 


Mw™ X 10^ 


2 


MuJrd X 10^ 


1 


0.992 


5.50 


0.76 


29.6 


2 


0.990 


5.50 


0.80 


29.2 


3 


0.989 


5.50 


0.85 


28.8 


4 


0.979 


5.40 


0.85 


28.4 


5 


0.976 


5.75 


0.85 


27.9 


6 


0.978 


5.90 


0.78 


27.4 


7 


0.986 


5.95 


0.85 


26.9 


8 


0.991 


6.20 


0.85 


26.4 


9 


0.992 


6.15 


0.84 


26.0 


10 


0.991 


6.25 


0.80 


25.7 


11 


0.993 


6.20 


0.85 


25.6 


12 


0.993 


6.20 


0.83 


25.7 


13 


0.993 


6.15 


0.81 


26.0 


14 


0.990 


6.25 


0.80 


26.5 


15 


0.983 


6.20 


0.80 


27.2 


16 


0.989 


5.95 


0.85 


27.9 


17 


0.985 


6.15 


0.80 


28.6 


18 


0.988 


5.82 


0.83 


29.1 


19 


0.984 


6.00 


0.78 


29.6 


20 


0.985 


5.75 


0.80 


29.9 


21 


0.984 


5.50 


0.83 


30.1 


22 


0.991 


5.50 


0.80 


30.1 


23 


0.992 


5.45 


0.77 


30.1 


24 


0.992 


5.47 


0.81 


29.9 


25 


0.992 


5.55 


0.84 


26.9 


26 


0.992 


6.05 


0.74 


29.0 


27 


0.993 


5.20 


0.85 


26.9 


28 


0.988 


5.30 


0.80 


23.6 



Table I: Values of the overlap integral given by eq.(10l between the analytical and the 28 test waveforms: the first 24 are 
numerically generated, with initial conditions Si/mj = (sina,0,cos a) and 82/^2 = (—0.6,0,0), with a — {k — 1) x 7r/12 for 
< fc < 24, the following 4 have spins parallel to the orbital angular momentum with components respectively {Siz, S2z)/m^ = 
(0, 0), (0.6, 0.6), (0.6, —0.6), (—0.6, —0.6). In the columns from third to fifth the values of ujrn and z maximizing the overlap and 
of ujrd are reported. 



ojm-Coeff. 


UJ,n fit (XIO^^) 


z-Coeff. 


z fit (xlO^^) 


ao 


55.500 


60 


84.00 


ai 


0.997 


bi 


-2.15 


13 


-2.032 


bs 


-4.42 


a4 


5.629 


bi 


-2.64 


15 


8.646 


bs 


-5.88 


ae 


-5.909 


fee 


-2.21 


aj 


1.801 






aj 


-14.059 






as 


15.483 






ag 


8.922 







Table II: Coefficients of eqs. (Ill and (131 as determined by comparison of the analytical waveforms with the 28 short test 
waveforms. 
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Figure 1: Quantitative comparison between the calibrated PhenSpin waveforms and the numerical relativity ones (top) and 
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